page blank not filmed 



38. Radial Pressure in the Solar Nebula as Affecting the 
Motions of Planetesimals 
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In the classical rotating Laplacian-type nebula, pressure gradients can develop 
radially to the protosun because of central radiation , particle ejection , and magnetic- 
field expansion or because of radial temperature or total gas density gradients. Except 
for the last two effects , the acting central acceleration for the gas is reduced from the 
gravitational value; the pressure gradient in the gas caused by temperature or density 
gradients may either add to or subtract from the gravitational acceleration , depending on 
the sense of the pressure gradient. Planetesimals in the nebula may thus experience 
tangential accelerations (4- or — ) with respect to the gas because of the differential 
radial accelerations acting on the particles and the gas. -4s a consequence } the plane- 
tesimals may spiral outward or inward with respect to the protosun. The present paper 
deals with growing planetesimals and a range of drag laws depending on the Reynolds 
number and on the ratio of particle size to mean free path. 

Particles spiral in the direction of positive pressure gradient, thus being concentrated 
toward toroidal concentrations of gas. The effect increases with decreasing rates of 
particle growth , i.e. f with increasing time scales of planet formation by accretion In the 
outer regions , where evidence suggests that comets were formed and Uranus and A eptune 
were so accumulated , the effect of the pressure gradient is to clear the forming comets 
from those regions. The large mass of Neptune may have developed because of this 
effect, perhaps Neptune’s solar distance ivas reduced from Bode’s “law” and perhaps 
no comet belt exists beyond Neptune. In the aster aid belt , on a slow time scale, the effect 
may have spiraled planetesimals toward Mars and Jupiter , thus contributing to the lack 
of planet formation in this region . 


BASIC PRINCIPLES 

L et us assume that during planetary 
formation in our solar system, gas and dust 
grains constitute a Laplacian-type nebular disk 
rotating about a central mass, roughly that of the 
Sun today. The grains or planetesimals are 
generally accreting in all size ranges from atomic 
and molecular to large bodies. For smaller bodies, 
the gas, presumed to be approximately a solar 
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mix, acts as a buffer against high velocities of 
encounter that might cause collisional destruc- 
tion. The motions arc fairly circular, controlled by 
a nearly inverse-square law of central force, while 
the gas and, at first, the grains are spread about 
the fundamental plane, held by the gas pressure 
against the surface gravity in the plane and the 
perpendicular component of the central force. 

In such a system the gas pressure varies from 
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point to point, depending on gravity and possible 
magnetohydrodynamic effects if a plasma is 
present. Radial variations in gas pressure affect 
the net central radial acceleration on the gas, 
causing the gas to deviate from the motion appro- 
priate to the local gravitational acceleration and 
therefore appropriate to the motions of the solid 
grains, as shown in figure 1. The grains thus meet a 
resisting medium and tend to move in toward the 
protosun or are accelerated radially outward, 
depending on the sense of the effective radial 
acceleration acting on the gas. 

Let A g be the deviation from the central gravita- 
tional acceleration acting on the gas, whether 
caused by spiraling magnetic fields or by a radial 
pressure gradient in neutral gas, dP/dr, where 
P is the total gas pressure and r is the radial 
distance to the center of the protosun. For a gas of 
density p, the radial correction to the central 
acceleration is, by the classical formula, 


Ag- 


1 dP 
^ p dr 


( 1 ) 


Note that in equation (1) the sign is reversed 
from that for stellar interiors because the gas in 
the nebula is held in radial equilibrium by its 



EFFECT OF GAS PRESSURE GRADIENT ON PARTICLE MOTION 


Figure 1 . — Effect of gas pressure gradient on particle 
motion. 


motion in the central gravitational field for 
dP/dr = 0. A positive pressure gradient outward 
from the center increases the effective central 
acceleration from that of gravity alone and causes 
the gas in equilibrium motion to rotate more 
rapidly. In that ease, the gas accelerates the grains 
in their near-circular motion and causes them to 
move radially outward. The opposite is true for 
dP/dr <0. 

To determine the quantitative motion of the 
grains, note that the near-circular velocity change 
AV in the gas under a central acceleration g 


follows from the velocity law 


y = 01/2,4/2 

(2) 

by differentiation, to the form 


r m Ag 

2g w 

(3) 


A spherical grain of radius s, density p a , and 
mass m will experience a drag ( — ) or accelerating 
( + ) force F by interaction with the gas at relative 
velocity v , where v<AV. The drag equation in a 
neutral gas is formally stated as 

Q 

F=» y xsV' 2 (4) 


where Cd is the dimensionless drag coefficient 
dependent on the Reynolds number R e if the mean 
free path L of the gas atoms or molecules is small 
compared to s . 

For gases the Reynolds number is defined by 



(5) 


where rj is the viscosity of the gas, given approxi- 
mately by i? = 3 ^pvL, in which v is the mean 
kinetic speed of the atoms or molecules. 

In case L<$, following the approximations by 


Probstein and Fassio (1969), 



C D = 24i?, _1 

for 

ff„<i 

(6a) 

C D = 24Rr m 

for 

KA.CIO" 

(6b) 

* 

II 

£ 

for 

R t > 10 3 

(6c) 


where equation (6a) represents Stokes’ law of drag 
and (6c) the Newtonian case, which holds fairly 
well for supersonic velocities (v>v). 

For the subsonic case when L>s, the drag 
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equation becomes generally complicated, but for 
v<Kv } the Epstein approximation has some 
validity: 

F(Epstcin) -Avps 2 v (7) 

where A— 47 r /3 if the accommodation coefficient 
of atoms or molecules on the moving sphere 
is unity. 

For our purposes, the various forms of the drag 
equation can be simplified by the introduction of 
the quantity stopping time t €i which is the ratio of 
the body's speed to the absolute value of the 
deceleration from drag. Equations (6a), (6b), 
(6c), and (7), respectively, can then be expressed 
as follows: 


R e < 1 



(6'a) 


1<F,<10 3 


2 3/5 pi$ 8/5 

9 ^/ 5^275 


Re> 10 3 



(6'c) 


Epstein (7') 

pv 

where for all, the drag law becomes simply 

( 8 ) 

TfX t e 

In the case of equations (6'b) and (6'c), t c has 
been changed from the e-folding decay time by a 
factor of less than 2 to produce the generality of 
equation (8) , which is correct for all cases, Stokes' 
Iaw T (eq. (6'a)) attains an error of only a factor of 
3 up to R e — 10 2 . 

The motion of the spherical grain or planetesimal 
in the solar nebula under influence of the gas 
“wind" is generally complicated. The present 
paper will deal in detail w T ith only two extreme 
and simple illustrative cases: (A) where the 
“drag” force F for r~AF is large so that the grain 
at velocity &V with respect to the gas would be 
stopped in a small fraction of the period of revolu- 
tion T Pi i.e., f««T P /2 tt; and (B) where the grain 
or planetesimal at relative velocity v — AF would 
not be stopped for a time comparable to or greater 
than the period of revolution, i.e., t e ^>T p /2r. 


In case A, where t e <&T p/2?r, the grain will 
largely follow^ the rotational motion of the gas and 
thus experience a central acceleration reduced or 
increased by Ag. Hence, it will execute a radial 
motion at velocity v with respect to the gas 
derived from F(v)c^m Ag. 

In case B, where t<y>T v /2w, the much larger 
grain or planetesimal will experience a tangential 
drag or an accelerating force at velocity ' 'AV 
and slowly spiral radially in or out with respect 
to the rotating system of gas according to the 
modified equations of motion in a central force. 


EQUATIONS OF MOTION 

Case A.— Grains carried along with the gas, 
t&Tp/ 2 tt. The relative velocity v of the grain in 
the gas is smaller than AV and very much smaller 
than the speed of sound. Thus, case A divides into 
only two subcases: case Al, when the mean free 
path L is larger than the particle radius s, requiring 
Epstein's law (eqs. (7), (7'), and (8)); and case 
A2, when L<s and R e is small, requiring Stokes' 
law r (eqs. (6'a) and (8)). 

Case Al — K<AF, v<Kv, L>s, and t e «T p /2r. 
For a grain (radius s, mass m, density p«) the 
terminal radial speed in the cloud under the 
acceleration Ag becomes from Epstein's law 

(eqs. (7') and (8)) 

dr A p*s /A v 

-j- —U Ag— — Ag (9) 

dt pv 

where Ag is given by eq. (1) for a radial pressure 
gradient in the gas or by the effect of a radial 
acceleration on the gas from some other source 
such as a magnetic field. 

Case A2— t«Al T , K<t), L<s f and t e «T p /2t. 
Here the terminal radial speed of the grain is 
derived from Stokes' tew (eqs. (6'a) and (8)), so 
that 


dr 2 p s s 2 Ag 

- =t,A g = — 

dt 


where r\ is the viscosity of the gas. 

Equation ( 10) is independent of the gas density 
so long as L<s and holds fairly well to = 10 and 
within a factor of 3 to Re — 10 2 , 

Case B. — Grains and planetcsimals meeting a 
resisting medium at v^A\ <3CF. Here, t e ^.T p and 
generally L<s, v<Kv t while R e may become 
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appreciable for large bodies, although in most cases 
of interest R e is small. Tisserand (1896) solved 
this two-body problem for small orbital eccen- 
tricity e and constant v-AV<KV. To the first 
order in e, de/ dt = 0, so that a small orbital eccen- 
tricity is retained either under resistance or 
acceleration for gas in nearly circular motion. 

For circular velocity V of mass m at distance r 
from the effective center of mass M, and a re- 
sisting or accelerating force F, the classical two- 
body solution takes the form 

1 dr 2F(AV) 2AV 1 Ag 

r dt mV t e V t e g K 

by equations (2), (3), (5) and (8). 

The appropriate expression for t e can be chosen 
from equation (6'a), (6'b), or (6'c), depending 
on the value of R t for v= AV (eqs. (5) and (4)). 
Since the spiraling rate is extremely small for 
large bodies with large values of R e , only the 
Stokes’ law (eq. (6'a)) equation will be presented 
explicitly, viz : 

dr 9 rrr Ag 

i‘w7 (12) 

The intermediate case between A and B, when 
t e ~T v / 2ir, has been solved by F, Franklin (private 
communication) for Epstein's and Stokes' laws, 
the drag being proportional to v. The maximum 
rate of orbital change is dr/dt = 0.5A7 at t e = 
0.8T p /2n . When t t — T P} dr/di equals 0.8 the value 
given by equation (11), approaching equation 
(11) as a limit for larger values of t e . For t e <T p / 30, 
equations (9) and (10) are correct within about 
20 percent and improve for smaller t e . The inter- 
mediate case thus gives values of dr/dt varying 
from 0.3AF at t e =T p / 30 to a maximum of 
0.5AF at t e =T p /8 and dropping to 0.3AF at 
t e = T v . Near but outside these limits, equations 
(9) and (10) or equation (11) gives fairly satis- 
factory values for dr/dt. 

EFFECTS OF GRAIN GROWTH ON 
RADIAL MOTION 

We wish to consider grains that are growing by 
the accretion of atoms and molecules from the 
nebula. For atoms or molecules that can “freeze” 
on the grain, constituting a fraction /» of the total 
density by weight, with an accommodation 


coefficient a iy a mean molecular velocity v % , and a 
mean diffusion coefficient D if a spherical grain of 
radius s and density p a grows at a rate 


ds _ otjftVip / ViS \ 1 
dt 4 p, \ + 4Z) t / 


(13) 


The second term becomes significant when the 
grain has reached such a size that for a velocity v 
through the gas small compared to grain growth 
is inhibited by gas diffusion of the appropriate 
atoms and molecules. Approximately, D,= lArji/p, 
where 77* is the viscosity appropriate to the “freez- 
ing” atoms and molecules. Hence, from the second 
term of equation (13), the grain growth rate slows 
when s'2>5.§7}j/vip, and is given by 

ds 1 Aoti f iVi 

Tt = — (14) 

dt p a s 

becoming independent of the general gas density 
and varying inversely as the radius of the grain. 

The transition from rapid grain growth (first 
term of eq. (13); e.g., Kuiper, 1951) to diffusion- 
limited grain growth (eq. (14)) begins when the 
mean free path becomes small compared to s, 
or roughly when the drag law at low grain veloci- 
ties changes from the Epstein law (eq. (9)) to 
Stokes' law (eq. (10) or (12)). Hence, for the 
Epstein case (Al) we keep only the first term in 
equation (13), so that a very small grain grows to 
radius Si in time h given by 


ti — 


4 p» ^ 

UifiViP 


(15) 


Again for the Epstein case, by equations (9) 
and (15), the change in radial distance from r 0 at 
t = 0 to n at h and $i becomes 


otifiViAg 2 Pa 2 Ag Sl 2 

1 0= T — tl = — 7^1 

8v on JiViVp 2 

applicable when v<^v, U<£T p /2t ) and L<s and 
when changes in p and g with r are neglected. 
Equation (16) applies well in the assumed low- 
density region of the solar nebula beyond Uranus, 
where cometary accretion is expected. Note that 
exhaustion of condensable gas slows the growth 
rates and thereby increases the total amount of 
spiraling. 

Within Jupiter's orbit, the relatively high gas 
density usually assumed and the consequent rapid 
growth rates carry the grain through the regimes 
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of Epstein's and Stokes* laws rather quickly . Thus, 
equation (12) can be used if applicable, and we 
will introduce the equations for large Reynolds 
numbers (t* from eq. (6'e) and dr/dt from eq. 
(11)) and larger bodies, even though the relative 
velocity AF (eq. (3)) may be small. On this basis, 
the slow growth from s 2 to s 3 wall take place by 
equation (14) in time t%—t 2 given by 


p a (sz 2 — s 2 2 ) 

U ~ h ~ 2 . 8 «,/.•* 


(17) 


while 


1 dr = pjAVy 1 
r dt 3 p s V $ 


(18) 


where the sign of dv/dt is that of AT , given bj 
equation (3), AF = 2^ 2 <r 1/2 Ag. 


APPLICATION IN THE 
“COMETARY REGION’' 

Let us arbitrarily define the “cometary region 7 * 
at r = 25 ATI, with a typical density in the plane 
of 10“ u 4 g cm" 3 , T = 55° K, P = 10~ 21 dyn cm“ 2 , 
giving a total areal density across the plane of 
10 2 3 g cm -2 (2X10" 5 solar mass per square AU) 
for a solar mixture by weight of gases (H, He, Nc, 
Ar) 0.9803, ices (C, N, 0 plus II) 0.0175, and 
Earthy materials (heavy elements) 0.0022. The 
corresponding viscosity will be 7? = 10 4 3 egs, 
fj= 10 4 85 cm s -1 for a mean molecular weight of 
2.34, and 10 4 41 cm s^ 1 for the ices plus Earthy 
material of mean molecular weight 18.4 and a 
mass fraction = 0.0197. 

The peak of the gas pressure and density will 
have been farther out, near r = 30 AU at Neptune's 
present solar distance, but the ordcr-of-mag- 
nitude calculations at r = 25 AU will illustrate the 
nature of the spiraling phenomena caused b^ a 
negative pressure gradient outward near the edge 
of the nebula. Suppose the pressure falls linearly 
with r to 0 at 50 AU. Then dP/dr =- 10“ 16 7 cgs, 
while Ag=10 -5,3 cm s~ 2 by equation (1) and 
AF=10 3 * 2 cm s _1 by equation (3) for a central 
solar mass. A forming icy grain of radius s and 
density 0.1 g cm -3 thus meets a resisting velocity 
of 16 m s -1 , small compared to molecular veloci- 
ties. The Reynolds number for equation (5) 
becomes R e = 10~ 3,6 s (cm), so that we are in the 
realm of Epstein’s and Stokes 7 laws of drag np to 


“eometesimals 77 of radius 100 m or larger at great 
solar distances. Correspondingly, the mean free 
path of the gas molecules is of the order of 1 m or 
more so that little error is made in applying 
Epstein's law up to perhaps s = 10 m, while to 
about the same limit the diffusion slowing of the 
accretion rate by the second term of equation 
(13) can be neglected. Hence, by equation (15), 
the time h for a tiny grain to grow to a radius $1 
becomes 

^ = 6.3sx(cm)yr (19) 

if we accept an accommodation coefficient «*= 1.0. 

Thus, our eometesimal will grow to a radius of 
10 m in some 10 4 yr. Its loss of solar distance from 
7*0 to 7*i, given by equation (16) , then becomes 

r 0 — 7*1 = 1 .2 X 10“ 5 si 2 AU (20) 

so that the eometesimal spirals in “nominally 
12 AU 77 by the time it has grown to a radius of 
10 m, neglecting significant changes in density 
and other quantities depending on solar distance. 

Note that the time varies inversely as the nebu- 
lar density p and that the amount of spiraling 
varies as p so that the inward spiraling for a 
constant Ag^dP/drp is increased for cometesi- 
mals of a given size at greater solar distances or at 
lower solar densities, even though the time scale 
increases as p~ l . Reduction in growth rates by any 
cause increases the total amount of spiraling for 
a given gas density and gradient. 

The sizable magnitude of the spiraling rate for 
eometesimals growing at the edge of the solar 
nebula may explain three facets of the present 
solar system: (a) the comparable masses of 
Neptune and Uranus, (b) the reduction of 
Neptune's solar distance from Bode 7 s “law, and 
(c) the still unobservable “comet belt beyond 
Neptune, expected by the writer (Whipple, 1964) . 

If we assume (Kuiper, 1951; Cameron, 1962; 
Whipple, 1964) that Uranus and Neptune are 
aggregates of eometesimals, then it is otherwise 
surprising that there should have been enough 
material at Neptune's distance to make a planet 
as large as Uranus. The spiraling effect, however, 
could bring in the eometesimals to a solar distance 
where the pressure gradient was smaller and pile 
up the cometary material from greater distances 
for accretion at Neptune, producing no sizable 
planet (Pluto?) beyond. Neptune's mean solar 
distance may also have been reduced. Thus, my 
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expectation of a comet belt beyond Neptune is 
perhaps unwarranted. In fact, the observations 
indicate that less than one Earth mass exists in a 
ring beyond Neptune to a solar distance of 50 AU 
(Hamid, Marsdcn, and Whipple, 1968). 

APPLICATION IN THE 

ASTEROIDAL REGION 

We may take a “typical” asteroidal region at 
t = 2.5 AU, with density in the plane of 10 -8 4 g cm~ 3 
or a pressure of 10~ 41 atm, T=550° K, P=I0 +1 9 
dyn cm 2 , giving a total areal density across the 
plane of 10 4 3 g cm -2 (2X1CL 3 solar mass per 
square AU), 77 = 10-^ s cgs ^ £j =10 5.35 cm ^ and 

Vi= 10 4 78 cm s- 1 for Earthy molecules of mean 
molecular weight 33.5 and a mass fraction /,= 
0.0022. The assumed temperature is that derived 
by Larimer and Anders (1967), supported by the 
further conclusions of Keays, Ganapathy, and 
Anders (1971) that in 13 chondrites at formation, 
7 7 =530 itgQO K and P=10"~ 4±2 atm. The atmos- 
pheric model is here calculated to allow for one 
central solar mass plus the additional surface 
pressure derived from the gravitational attraction 
of the gas on itself. 

Suppose, for illustration, that the gas pressure 
doubles in 1 AU toward the Sun so that dP/dr = 

— 10~ 113 cgs, A<7= — 10~ 2 - 9 cm s~ 2 by equation (1), 
and AF = 10 4 - 1 cm s _I by equation (3). In the gas 
the mean free path is 0.2 cm and R e = 0.64s. This 
puts us squarely in the Newtonian range of drag 
(eq. (6'c)) and the range of slow grain growth 
(eq. (15)) for sizable planetesimals of $>10 m. 
Growth to such sizes for p* = 3 occurs rapidly, 
~ 1.4 cm yr' 1 , for a short time and then slows. As 
we have seen, the radial spiraling rate for a change 
approximately fifty fold in radius by Stokes’ 
law averages about 0.3 AV or dr/dt= ICC 2 * 1 AU 
yr -1 . Thus, the reduction in r is a fraction of 1 AU 
while the planetesimal is growing to meter dimen- 
sions. 

For larger bodies, $>10 m, we can apply 
equation (18) (Newtonian drag and the classical 
resisting medium) to find the spiraling distance 
Ar in time At for constant s (cm) 

Ar 1.22 

7 =~— M(yr) ( 21 ) 

Appreciable change of solar distance occurs 


when At In years becomes comparable to the radius 
of the planetesimal measured in centimeters and 
varies directly as the nebular density (eq. (18)). 
The growth rate, however, is really quite un- 
certain. The assumed value of AF is still small 
(less than 1 percent) compared to the orbital 
velocity but is some 7 percent of the mean 
molecular velocity. Turbulent motion of the gas 
might easily produce random relative velocities 
of the order of 1 km s” 1 so that the faster growth 
rate might be applicable (first term of eq. (13)). 
In that case, the change of r for a further growth 
to kilometers might not be significant. The slow 
growth rate, however, would make a significant 
change in r quite likely. 

We may conclude, therefore, that effects of 
pressure gradients in the asteroid belt could be 
very important in shifting planetesimals toward 
Mars or Jupiter, provided the growth rates of the 
asteroids are not too rapid. For the basic develop- 
ment of the asteroids on a short time scale of 10 3 
to 10 4 yr, the effect would be minor. However, 
should the nebular density be 10~ 5 atm or less in 
the asteroid region and the growth time 10 6 to 
10 8 yr, pressure gradients could well have deci- 
mated the asteroid region by spiraling the plane- 
tesimals toward Mars or Jupiter, depending on 
the distribution of gas and the location of original 
grain growth. Possibly this effect contributed to 
the failure of an Earth-sized planet to develop 
between Mars and Jupiter. 

GENERAL REMARKS 

My attention was first directed to systematic 
interactions between gas motions and growing 
planetesimals in a postulated solar nebula by 
Hoyle’s (1960) theory involving the expansion of 
a centrally condensed nebula by the outward 
force of spiraling magnetic fields from the Sun. In 
this process, the plasma w-ould be pressed outward 
from the Sun, carrying with it the neutral gas and, 
according to Hoyle, also the planetesimals. In fact 
(Whipple, 1964) , the reduced effective gravity on 
the gas will, following the arguments of this paper, 
cause the planetesimals to spiral inward toward 
the Sun rather than outw f ard. 

Cameron (1969) mentions the effect of a radial 
pressure gradient in producing an inw ard spiraling 
of planetesimals but docs not discuss the alterna- 
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tive possibility, viz, outward spiraling, should a 
positive outward pressure gradient occur. Gen- 
erally, if a toroid of higher density occurs in the 
solar nebula, the growing planetesimals are drawn 
toward it from the inside as well as from the 
outside, increasing the growth rate of an accreting 
planet. The importance of the effect, as we have 
seen above, depends mostly upon the nebular 
density and the time rates of chemical condensa- 
tion. For a short time scale (~1(P to 10 4 yr) such 
as Cameron envisages, the pressure-gradient 
effect would not be important. 

For long time scales of planetary formation, 
10 6 to 10 s yr, as are frequently postulated, the 
pressure-gradient effect could be highly signifi- 
cant. It may well have drawn in the cometary 
material from far beyond Neptune, adding 


materially to Neptune's mass and reducing the 
solar distance from Bode's “law." It may have 
greatly reduced the comet production in these 
regions of space so that no significant comet belt 
exists beyond Neptune. Conceivably, the pressure- 
gradient effect may have assisted in the building 
of Mars and Jupiter, at the expense of the asteroid 
belt. A lower surface gas density between the 
regions of Mars and Jupiter might not, in itself, 
have prevented the accretion of a sizable terrestrial 
planet between Mars and Jupiter without the 
action of the pressure-gradient effect. 
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